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Present applications of the dispersive-optical-model analysis are restricted by the use of a local 
but energy-dependent version of the generalized Hartree-Fock potential. This restriction is lifted by 
the introduction of a corresponding nonlocal potential without explicit energy dependence. Such 
a strategy allows for a complete determination of the nucleon propagator below the Fermi energy 
with access to the expectation value of one-body operators (like the charge density), the one-body 
density matrix with associated natural orbits, and complete spectral functions for removal strength. 
The present formulation of the dispersive optical model (DOM) therefore allows the use of elastic 
electron-scattering data in determining its parameters. Application to 40 Ca demonstrates that a fit 
to the charge radius leads to too much charge near the origin using the conventional assumptions 
of the functional form of the DOM. A corresponding incomplete description of high-momentum 
components is identified, suggesting that the DOM formulation must be extended in the future to 
accommodate such correlations properly. Unlike the local version, the present nonlocal DOM limits 
the location of the deeply-bound hole states to energies that are consistent with (e,e'p) and (p,2p) 
data. 

PACS numbers: 21.10.Pc,24.10.Ht,11.55.Fv 



I. INTRODUCTION 

An important link between the description of nuclear 
reactions and nuclear structure was proposed in Ref. [l[ , 
establishing a workable connection between the optical- 
model and the shell-model potential. The proposed 
implementation is now known as the dispersive opti- 
cal model (DOM) and has been extensively reviewed in 
Ref. [2[. Using empirical knowledge of standard opti- 
cal potentials Mahaux and Sartor proposed to employ 
the experimentally well-constrained imaginary parts of 
the optical potential as the critical building blocks to 
determine, through a subtracted dispersion relation, the 
corresponding dynamic real part without requiring addi- 
tional parameters. Since the subtraction point is usually 
chosen to be the Fermi energy, the remaining static and 
real Hartree-Fock-like potential can be linked to empir- 
ical information from mean-field theories. Conventional 
Woods-Saxon form factors have been employed includ- 
ing standard forms for volume and surface contributions. 
By assuming a similar energy dependence of the imagi- 
nary part of the potential above and below the Fermi en- 
ergy, it is then possible to derive and successfully predict 
properties of bound nucleons, including removal energies, 
overlap functions, spectroscopic factors, and widths [2j. 

The success in reproducing spectroscopic factors de- 
rived from the (e 7 e'p) reaction [3j relies partly on the 
inclusion of an energy asymmetry far from the Fermi en- 
ergy for the imaginary volume contribution associated 
with the larger phase space of particle-like as compared 
to hole-like states. This asymmetry can be inferred from 
nuclear-matter calculations of the imaginary part of the 
self-energy A number of studies have made use of 



the framework of the DOM (34l2j. Some steps towards 
a global version of the DOM have been recently reported 
in (T3 - [lpl. A different perspective was recently devel- 
oped in [l7l Il8j where the DOM was exploited to extract 
the nucleon asymmetry dependence from fitting 40 Ca and 
48 Ca data. This allows both interpolation and extrapo- 
lation to other nuclei and the prediction of experimental 
data. Such data-driven extrapolations may play an im- 
portant role in predicting the properties of nuclei towards 
the drip lines. 

Additional benefits may result when DOM ingredients 
are used in the description of transfer reactions using the 
adiabatic distorted-wave approximation which employs 
proton and neutron optical potentials for the description 
of the relevant deuteron scattering wave function [l9|, H(| ■ 
The spectroscopic factors deduced from transfer reac- 
tions continue to exhibit a strong dependence on which 
optical potential is employed, as was recently shown for 
a number of Ar isotopes [2l[ . In addition to this uncer- 
tainty, there is a substantial discrepancy between the ex- 
tracted single-particle (sp) properties from transfer and 
heavy- ion knockout reactions [22, 23] . Without a means 
to unambiguously extract such sp properties, the success 
of rare-isotope facilities will be severely hampered. A 
well-constrained complete optical potential in the sense 
of the self-energy for the Dyson equation for nucleons will 
therefore be an important ingredient in obtaining unique 
and undisputed information [24J linking structure and re- 
action data in a unified manner. 

The DOM can easily be applied to elastic-scattering 
data, since only phase shifts are required to describe dif- 
ferential cross sections and polarization data. To be use- 
ful for transfer reactions, knowledge of the interior wave 
functions of protons and neutrons at positive energy is 
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required. The current implementation of the DOM em- 
ploys a real component, a Hartree-Fock-likc term, that is 
in principle nonlocal but is replaced by an equivalent lo- 
cal potential with an energy dependence mainly governed 
by the so-called A;-mass [2] . Since this energy dependence 
does not result from a dispersion integral, it leads to a 
distortion of the normalization that can be approximately 
fixed [25l - l27l |. Since specific assumptions about a gaus- 
sian form of the nonlocality are made in addition to the 
approximate nature of the local approximation, it is quite 
important to explore the actual form of nonlocal poten- 
tials that are based on microscopic calculations of the 
nucleon self-energy like in the Faddeev random phase ap- 
proximation (FRPA) [28-30]. Another reason to explore 
the inclusion of a nonlocal Hartree-Fock-like potential is 
to allow additional data to be included in the fitting pro- 
cedure. It is the goal of the present paper to explore the 
inclusion of a nonlocal Hartree-Fock-like potential in the 
DOM in order to describe a larger set of data in partic- 
ular those pertaining to properties of nucleons below the 
Fermi energy. 

The normalization distortion is particularly significant 
below the Fermi energy, where it leads to difficulties in 
determining the sp strength distribution A nonlocal 
DOM potential can avoid these problems and provide 
a properly normalized solution of the Dyson equation. 
The resulting propagator below the Fermi energy then 
provides access to spectral functions, the one-body den- 
sity matrix, and all one-body expectation values in the 
ground state. In turn, elastic electron-scattering data 
that yield the nuclear charge distribution can be used 
to constrain the DOM potentials. The sp strength dis- 
tributions obtained from (p,2p) [HI, HH and (e,e'p) Q 
reactions provide further constraints. A nonlocal DOM 
potential also avoids the pitfall of the linear energy de- 
pendence of the equivalent local potential which provides 
too much binding for the most deeply bound levels. 

Experimental evidence for the presence of high- 
momentum components in the nuclear ground state has 
been provided by the (e 7 e'p) reaction 33] in reason- 
able agreement with microscopic calculations for light nu- 
clei [Mill. While the number of high-momentum pro- 
tons in light nuclei represents a modest 10% of the total, 
their presence does confirm the basic tenet of most real- 
istic nucleon-nucleon interactions that contain a sizable 
repulsion at short relative distances. It is unclear whether 
current DOM implementations generate any sizable pres- 
ence of high- momentum nucleons. It is therefore useful to 
explore the spectral distribution in momentum space of 
the DOM propagator to assess its high-momentum con- 
tent, and, if found lacking, provide an incentive to con- 
struct DOM self-energies that represent the experimen- 
tal findings (33|. While spectroscopic factors obtained 
from the analysis of the (e,e'p) reactions are already em- 
ployed in current DOM fits [IT], Ell], a nonlocal treat- 
ment of the DOM potential opens the possibility to use 
these cross-section data directly in the fitting procedure, 
since both the overlap function of the removed proton as 



well as the scattering wave function of the outgoing pro- 
ton can be described by the corresponding DOM wave 
functions. Data from the (e,e'p) reaction therefore can 
provide additional constraints and provide further con- 
firmation of the interpretation of these data as well as 
reducing the uncertainty in the quoted absolute spectro- 
scopic factors [24j |. 

The purpose of this work is to clarify the inclusion 
of a nonlocal Hartree-Fock-like potential and generate 
the resulting solution of the Dyson equation below the 
Fermi energy, while keeping the earlier DOM results ob- 
tained with the local equivalent potential intact or im- 
proving upon them. We will continue to make the sim- 
plest possible assumptions about the form of the nonlo- 
cality and thereby not introduce any additional param- 
eters. In Sec. |TT]we will clarify the required strategy to 
solve the Dyson equation in coordinate space while pre- 
serving the ingredients obtained from earlier fits to 40 Ca 
data [Tt], [HI . The results illustrating the complete solu- 
tion of the Dyson equation below the Fermi energy are 
illustrated and discussed in Sec. IIIII A summary and 
conclusions are presented in Sec. |TV] 



II. DYSON EQUATION AND DOM 
SELF-ENERGY WITH NONLOCAL 
POTENTIALS 

A. Green's function ingredients 

We start with a brief summary of relevant results from 
the Green's function formulation of the many-body prob- 
lem (27| ■ The nucleon propagator with respect to the 
A-body ground state is given by 
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where complete sets of states in the A ± 1 systems are 
inserted and the sp basis with good radial position, or- 
bital angular momentum (parity) and total angular mo- 
mentum is chosen while suppressing the projection of the 
total angular momentum and the isospin quantum num- 
bers. The continuum solutions in the Ail systems are 
also implied in the completeness relations. The numer- 
ators of the particle and hole components of the propa- 
gator represent the products of overlap functions associ- 
ated with adding or removing a nucleon from the A-body 
ground state. The standard development of Green's func- 
tion theory relates the nucleon propagator to the self- 
energy yielding the Dyson equation in the following form 



G ej (r,r';E) = G { " ) (r,r';E) 



(2) 



+ Jdr r 2 Jdf> r' 2 G i ff(r,r;E)j: ej (r,r'- 1 E)G ej (r',r'-E). 

For the present discussion the noninteracting propaga- 
tor involves only kinetic energy contributions. The nu- 
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cleon self-energy contains all linked diagrammatic contri- 
butions that are irreducible with respect to propagation 
represented by G^°> . All contributions to the propagator 
are then generated by the Dyson equation itself. The so- 
lution of the Dyson equation generates all discrete poles 
corresponding to bound A± 1 states explicitly given by 
Eq. JD that can be reached by adding or removing a 
particle with quantum numbers rtj. The hole spectral 
function is obtained from 



Sij(r\E) = -Im Gtj(r,r;E) 



(3) 



for energies in the continuum. The total spectral strength 
at E for a given £j combination, 
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yields the spectroscopic strength per unit of energy. For 
discrete energies as well as all continuum ones, overlap 
functions for the addition or removal of a particle are 
generated as well. The connection between the nucleon 
propagator and elastic-scattering data can therefore be 
made explicit by identifying the nucleon elastic-scattering 
T- matrix with the reducible self-energy obtained by iter- 
ating the irreducible one to all orders with [23, HE - 

m. 

For discrete states in the A — 1 system one can show 
that the overlap function obeys a Schrodinger-like equa- 
tion [13] • Introducing the notation 
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for the overlap function for the removal of a nucleon at r 
with discrete quantum numbers I and j, one finds 
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and in coordinate space the radial momentum operator 
is given by p r = —ih(-g- + i). Discrete solutions to 
Eq. ([6]) exist in the domain where the self-energy has 
no imaginary part and these are normalized by utilizing 
the inhomogeneous term in the Dyson equation. For an 
eigenstate of the Schrodinger-like equation [Eq. ©], the 
so-called quasihole state labeled by a q h, the correspond- 
ing normalization or spectroscopic factor is given by [27j 



S£ = 1 - 



dT,ij(a q h,a q h;E) 



dE 
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which is the discrete equivalent of Eq. Q . Discrete solu- 
tions in the domain where the self-energy has no imagi- 
nary part can therefore be obtained by expressing Eq. ([5]) 



on a grid in coordinate space and performing the corre- 
sponding matrix diagonalization. Likewise, the solution 
of the Dyson equation [Eq. (jSJ] for continuum energies 
in the domain below the Fermi energy, can be formu- 
lated as a complex matrix inversion in coordinate space. 
This is advantageous in the case of a nonlocal self-energy 
representative of all microscopic approximations that in- 
clude at least the Hartree-Fock approximation. Below 
the Fermi energy for the removal of a particle 
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the corresponding discretization is limited by the size of 
the nucleus as can be inferred from the removal ampli- 
tude given in Eq. ([5]). which demonstrates that only co- 
ordinates inside the nucleus need to be considered. Such 
a finite interval therefore presents no numerical difficulty. 



B. Link with the DOM 

While microscopic calculations of the nucleon self- 
energy have made substantial progress in recent 
years [28U301 . l39j . accurate representations of elastic- 
scattering data in particular have not yet gone beyond 
phenomenological representations of the self-energy in 
terms of standard or dispersive optical potentials. A clear 
link with the microscopic self-energy is provided by the 
DOM strategy [TJ, HJ. It employs the dispersion relation 
between the real and imaginary part of the microscopic 
self-energy given by 

, , , . , i\ f°°dE'ImZei(r,r';E') 
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where V represents the principal value. The static con- 
tribution arises from the correlated Hartree-Fock term 
and the dynamic parts start and end at corresponding 
thresholds in the A±l systems that have a larger sep- 
aration than the corresponding difference between the 
Fermi energies for addition (ej = Eq +1 — Eq) and re- 
moval (sp — Eq— Eq -1 ) of a particle. The latter feature 
is particular to a finite system and generates possibly sev- 
eral discrete quasiparticle and hole-like solutions of the 
Dyson equation in Eq. ([5]) in the domain where the imag- 
inary part of the self-energy vanishes. 

The standard definition of the self-energy requires that 
its imaginary part is negative, at least on the diagonal, 
in the domain that represents the coupling to excitations 
in the A + 1 system, while it is positive for the coupling 
to A — 1 excitations. This translates into an absorptive 
potential for elastic scattering at positive energy, where 
the imaginary part is responsible for the loss of flux in the 
elastic channel. It is convenient to introduce the average 
Fermi energy 
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Subtracting Eq. ([10]) calculated at this energy, from 
Eq. (fTU)) generates the so-called subtracted dispersion re- 
lation 
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The beauty of this representation was recognized by Ma- 
haux and Sartor [l|, [2| since it allows for a link with em- 
pirical information both at the level of the real part of 
the nonlocal self-energy at the Fermi energy (probed by a 
multitude of Hartree-Fock calculations) and also through 
empirical knowledge of the imaginary part of the optical 
potential (constrained by experimental data) that con- 
sequently yields a dynamic contribution to the real part 
by means of Eq. (|12p . This procedure requires further 
assumptions since detailed knowledge of the imaginary 
part of the self-energy below the Fermi energy has only 
become available with electron-induced proton knockout 
reactions Q . Since the empirical knowledge has relied on 
local representations of the imaginary part of the optical 
potential, it is natural to make a similar assumption for 
the DOM version. In addition, a separation in terms of 
surface (low energy) and volume (higher energy) absorp- 
tion has been incorporated in accordance with standard 
practice. 

Most implementations of the DOM simply assume that 
near the Fermi energy there is a similar behavior above 
and below the Fermi energy for the imaginary part of 
the surface contribution to the self-energy although this 
assumption requires more stringent tests. Farther away 
from the Fermi energy the phase space asymmetry char- 
acterized by the density of two-particle-one-hole (2plh) 
and two-hole-one-particle (2hlp) states leads to an as- 
sumed energy asymmetry for the volume contribution 
that is also consistent with expectations from nuclear- 
matter results. Standard forms and fitted parameters for 
a recent implementation of the DOM for Ca isotopes can 
be found in Ref. [18[ . 

A clear success of the DOM approach is the accurate 
prediction of bound sp properties including radii, spec- 
troscopic factors, and the compression of sp levels near 
the Fermi energy that originates from the dispersive con- 
tribution in Eq. (|12|) |2| and cannot be obtained with 
standard Hartree-Fock approaches. The latter feature 
provides an explanation for the observed nuclear level 
density parameter a which is large compared to sp esti- 
mates [40j . 

The remaining ingredient in the DOM procedure is the 
nonlocal self-energy at the Fermi energy represented by 
ReT,£j(r, r'\ Sf)- It is more convenient to discuss the spin- 
independent part of this term in a basis with vectors in 
coordinate space using the notation T,HF(f,r') and em- 
ploying the HF label that was introduced by Mahaux and 
Sartor Q even though this term is not a true Hartree- 



Fock contribution as the derivation of Eq. (IT2|) clarifies. 
The usual treatment of E#i?(r,r') is to assume that it 
can be rep laced by a local but energy-dependent poten- 
tial 0, |25| - |27| . The corresponding form then can be writ- 
ten as 



X HF (r,r') V HF (r,E)S{r-r' 



(13) 



where 



V H F(r,E) = V HF (E)f(r,r HF ,a HF ) (14) 
containing the Woods-Saxon form factor 
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The factorized linear energy dependence can be 
parametrized by the corresponding effective mass below 
the Fermi energy and can be represented by 



V HF (E) = Vhf{s f ) 



1 - 



L H F 



(E-e F ), (16) 



which can be combined with the Woods-Saxon form fac- 
tor to generate m* HF {r). This version is inspired by the 
Skyrme implementation of the HF potential [2j. More 
generally, one may identify this effective mass with an 
energy-dependent version of the effective mass m* (r; E) 
that governs the nonlocality of the self-energy and is 
sometimes referred to as the fc-mass. It was shown in 
Ref. [4l[ that this effective mass is critical to reconcile 
the phenomenological (local) imaginary part of the opti- 
cal potential with the microscopic one |27j and to explain 
the observed nucleon mean free path. For finite nuclei, 
this implies that the DOM version of its local imaginary 
part W is related to the self-energy by 

W (r; E) = m *^ ;E \ mE(r; E). (17) 
m 

This suggests that the use of a nonlocal HF self-energy 
in the DOM framework has to be accompanied by en- 
hancing the imaginary part with a corresponding factor 
m/m*(r; E). Results discussed later indeed corroborate 
the necessity of including this factor e.g. to obtain identi- 
cal spectroscopic factors as those in Ref. [l|| . It is there- 
fore possible to employ the same parameters as in the fit 
of [18| and only replace the energy-dependent local equiv- 
alent HF potential by a suitable energy-independent non- 
local one. We have chosen the standard form introduced 
in Ref. [H| to represent 

£ HF (r,r') - V NL f(h\r + r'\,r NL ,a NL )H(\r-r'\), 

where the degree of nonlocality is expressed by a gaussian 
governed by the parameter (3 
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This nonlocal form requires four parameters 
(YnLi r NLj O'NLj and /?), which is the same number 
required to represent VHF(r,E) in Eq. (fT4|) . We reit- 
erate that this nonlocal representation is essential in 
obtaining properly normalized spectral functions and 
spectroscopic factors. In the following we will discuss 
results for 40 Ca with this nonlocal version of the DOM 
with emphasis on energies below the Fermi energy. We 
note that all DOM parameters that were obtained in 
Ref. [HI have been kept while only the local HF potential 
with its spurious energy dependence has been replaced 
by Eq. ([TBI and the application of the dispersion relation 
Eq. (fT2|) has been modified according to Eq. (fT7|). 



III. RESULTS 
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As discussed in the previous section, only four parame- 
ters are needed to introduce the nonlocal HF contribution 
to the DOM. As in the usual DOM fit, the location of 
the main fragments of the 0d 3 / 2 and Is 1/2 valence hole 
levels was used to constrain the parameters of the non- 
local HF potential. Since the complete one-body density 
matrix can be obtained with a nonlocal HF potential, 
it was also possible to constrain the parameters by the 
mean square radius of the charge distribution that is well 
known experimentally [42j | . An additional problem that 
can be cured by the nonlocal version of the HF poten- 
tial is associated with the linear energy dependence of 
the local version as shown in Eq. (fl6|) . Typical DOM 
fits generate rather deeply bound Osi/2 states, often well 
below the peaks seen in (e,e'p) and {p,2p) experiments. 
With a nonlocal potential it is possible to use the peak of 
the deeply bound 0s 1 / 2 state as an additional constraint 
and avoid the problem. The resulting parameters are col- 
lected in Table U including those for the local potential. 
All other parameters and the detailed shapes chosen for 
the imaginary parts of the DOM potential for 40 Ca can 
be found in Ref. [T3|. When adjusting the parameters of 
the nonlocal potential it was found that it was possible to 
incorporate the constraint of the mean square radius of 
the charge distribution while generating quasihole frag- 
ments at energies that are at least as good as the original 
fit. It may be useful in the future to explore other forms 
for nonlocal potentials, especially when microscopic self- 
energies obtained with the FRPA method [28M30I [H| are 



TABLE I. Parameters for the local energy-dependent Woods- 
Saxon potential and the nonlocal version with gaussian non- 
locality for 40 Ca. 





local 


nonlocal 


Depth [MeV] 


-56.5 


-92.3 


Radius [fm] 


1.19 


1.05 


Difhiscness [fm] 


0.70 


0.70 


m* HF /m 


0.57 




Nonlocality [fm] 




0.91 



FIG. 1. (Color online) Comparison of proton S1/2 spectral 
strength with the nonlocal (solid) and local potential with 
the spurious energy dependence (dashed). The nonlocal result 
does not yet contain the correction given in Eq. (|17[) . Note 
that the dashed curve even exceeds the number of mean-field 
S1/2 particles by more than 50% illustrating the incorrect nor- 
malization when the local energy-dependent potential is used 
in the Dyson equation without proper corrections. 

analyzed. It is possible to keep the same value for the dif- 
fuseness parameter as for the local potential. Depth and 
radius parameters are however clearly very different and 
must combine with the nonlocality parameter /3 to gen- 
erate similar results for the energies of the valence hole 
states. It was found that the charge radius of the nucleus 
provides a significant constraint and should therefore be 
used in future applications of the DOM method. 

After projecting the nonlocal potential on to states 
with good orbital angular momentum, it is possible to 
perform the complex matrix inversion in coordinate space 
to solve Eq. ([2]) . We note that the imaginary part of the 
DOM potential of Refs. [l7|,[ll] ends at the Fermi energy 
[see Eq. (1TT1) ] so this procedure generates sharply peaked 
features for valence hole states just below the Fermi en- 
ergy. To illustrate the influence of the spurious energy 
dependence included for a local potential, we display in 
Fig. Q] the proton S1/2 spectral strength [see Eq. with 
a nonlocal potential (solid) and the local version (dashed) 
from the original fit. The nonlocal potential yields the 
same valence quasihole energy but the correction implied 
by Eq. (|T7)) was not yet applied. Spectral functions ob- 
tained by solving the Dyson equation with the spurious 
energy dependence associated with the local HF poten- 
tial typically overestimate the mean-field occupation by 
more than 50% leading to more than 30 protons for 40 Ca. 

While general features of the solid curve in Fig. [1] for 
S1/2 strength distribution appear in order, like the wide 
peak for the lowest orbit and a sharp well-localized peak 
for the one near the Fermi energy, it should be noted 
that the strength in the peak near the Fermi energy con- 
tains a spectroscopic strength of 0.78 not in agreement 



6 




I 1 1 I 

-90 -60 
E [MeV] 



E [MeV] 



FIG. 2. (Color online) Spectral functions for all mostly occupied £j combinations in 40 Ca 
together with the / 7 / 2 result. Orbits with j — £ + i are dashed to distinguish them from the ones with j =£ — ■ 
results exhibit similar peak locations and widths as observed in (p,2p) pjll. 1321] and (e, e'p) experiments [43 
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with the DOM results of Refs. [U El]. It is only when 
the remainder of the DOM potential is multiplied by the 
effective mass correction of Eq. (fTT)) that we obtain a 
spectroscopic factor of 0.65 for the lsi/2 state in bet- 
ter agreement with Ref. [l8j]. The correct strength dis- 
tributions for the most relevant proton £j combinations 
including fa/ 2 are displayed in Fig. [2] All peaks corre- 
spond to the orbits that are expected to be fully occupied 
in the mean field. It is important to realize however, that 
at each energy the total strength according to Eq. (j4]) is 
obtained for a given £j combination. For a given (reason- 
able) mean field potential, several of the corresponding 
orbits may exhibit finite amounts of strength at one en- 
ergy, including those that are not occupied in that mean 
field. The converse is illustrated by the fa/2 strength 
distribution, since this orbit is completely empty in a 
mean-field picture. The strength exhibited appears on 
account of the presence of the imaginary part of the self- 
energy below the Fermi energy which allows some finite 
amount of fa/ 2 strength to appear there. In addition to 
this feature, mostly occupied £j combinations exhibit a 
broadening of the strength with a width that represents 



the local mixing with more complicated states like 2hlp, 
etc. We note that it is straightforward to place the first 
peak of the s\/ 2 distribution in accordance with experi- 
mental observations from (p,2p) [3l|, [HJ and (e,e'p) ex- 
periments [43j]. This is more difficult when the local HF 
potential is employed due to its linear energy dependence 
which lowers the well with decreasing energy. When in- 
tegrating the total strength shown in Fig. [5] for all orbits 
except the fa/2 and multiplying with the corresponding 
degeneracy factor of 2j+l, the summed strength is 19.48. 
While this may appear reasonable, it should be kept in 
mind that the assumed state independence of the DOM 
potential (apart from spin-orbit) and the ^-dependence 
of the nonlocal HF potential (on account of its angu- 
lar dependence implied by Eq. (ITS)) ), implies that some 
strength will also be generated for higher £- values leading 
to an overestimate of the total proton number. Indeed, 
when the cut-off is placed at £ = 3, i.e. the fa/ 2 and f 5 / 2 
contributions are also included, the total proton number 
becomes 21.43. This suggests that in future DOM work 
the total proton (neutron) number should be used as a 
further constraint on the potentials. The possibility of 
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TABLE II. Quasihole energies for proton orbits in Ca for 
the local and nonlocal DOM implementation compared with 
experiment. 







Energy [MeV] 






orbit 


local 


nonlocal 


peak 


experiment 


OSl/2 


-57.3 


-47.4 


-46.7 


~-47 


OP3/2 


-35.1 


-31.4 


-31.1 


~-30 


Opi/2 


-30.3 


-26.7 


-26.4 


~-30 


0d 5/2 


-13.5 


-13.8 


-13.5 


-13.5 


1*1/2 


-9.5 


-9.8 


-9.8 


-10.8 


0d 3 /2 


-8.3 


-8.3 


-8.3 


-8.3 



TABLE III. Spectroscopic factors for proton orbits in Ca 
for the local and nonlocal DOM implementation. 



orbit local nonlocal 

OsL/2 LIT 0.98 

0p 3/2 0.94 0.93 

Opi/a 0.95 0.94 

0d 5/2 0.83 0.86 

ls 1/2 0.67 0.65 

0d 3/2 0.65 0.64 



including some state dependence may also be explored, 
in particular by rely in g on microscopic input from FRPA 
calculations [2i|3fl |39|| . 

Before discussing new results not available with the 
standard DOM implementation, we first compare sev- 
eral quasiparticle properties in the two approaches. In 
Table [TT] we compare quasihole energies obtained with 
the local and nonlocal DOM with experimental data (for 
deeply-bound orbits we use Ref. [32|). The numbers 
quoted in the following tables for quasiparticle properties 
using the local version of the DOM differ slightly from 
the ones generated in Ref. [l8|], since a small error in 
the calculation of the dispersive volume contribution has 
been corrected. The column labeled "local" reports the 
solutions of the eigenvalue equation for the local DOM 
potential without the imaginary part. This includes a 
self-consistency procedure since the potential is energy 
dependent, i.e. the chosen input energy has to coincide 
with the obtained eigenvalue. Such a calculation for the 
nonlocal DOM is reported in the column labeled nonlo- 
cal. The imaginary part is included in the column la- 
beled "peak" which identifies the location of the peak of 
the spectral function for each orbit. As expected, there 
is little difference between the latter two approaches, es- 
pecially close to the Fermi energy. The largest difference 
between the local and nonlocal approach occurs for the 
lowest S1/2 orbit. As discussed earlier, the nonlocal po- 
tential is better able to constrain the peak of the spec- 
tral S1/2 strength to the correct value. Overall agreement 
for the peak location appears quite satisfactory for the 
nonlocal DOM although it may be necessary to consider 
some state dependence if a better fit for Is i/ 2 quasihole 
energy is deemed appropriate. The results for neutrons 
are naturally not very different apart from the Coulomb 
shift to those for protons on account of isospin symmetry 
and are therefore not reported. We only note that the 
lowest neutron S\j 2 peak occurs at -56 MeV in the non- 
local DOM while in the local counterpart it is found at 
-67 MeV, confirming the discussion for the corresponding 
proton level that the local DOM tends to bind this orbit 
too deeply. 

In Table [TTT1 the spectroscopic factors are listed for the 
same orbits as in Table [II] These results are obtained 
for the local DOM by using the approximate expression 



for the spectroscopic factor reviewed in Ref. [2| and given 
explicitly in Eq. (12) of Ref. [l8j]. This expression does 
not guarantee that the resulting spectroscopic factor is 
less than 1 (as it should be), which is illustrated by the 
outcome for the 0si/ 2 orbit. For the nonlocal DOM, 
Eq. ([5]) has been used where the derivative is taken at 
the eigenvalue obtained from Eq. © with neglect of the 
imaginary part of the potential. This procedure is also 
not appropriate in the domain where the imaginary part 
becomes substantial and is already suspect for the d 5 / 2 
orbit. When the imaginary part is neglected, it is pos- 
sible that the total real dispersive correction has a pos- 
itive derivative at the energy corresponding to the self- 
consistent eigenvalue even in the nonlocal case, leading to 
an unphysical spectroscopic factor. Already the strength 
content of the peak for the d 5 / 2 orbit in Fig. [2] is more in 
line with the spectroscopic factors quoted for the lsi/2 
and 0d 3 / 2 orbits and therefore substantially smaller than 
the 0.86 listed in Table IIIII Only for the latter two or- 
bits is the neglect of the imaginary part of the poten- 
tial unimportant, since the content of the sharp peaks in 
Fig. [2] coincides with the spectroscopic factors given in 
Table IIIII We also note that there is reasonable agree- 
ment with the local and nonlocal DOM results for these 
levels. We therefore conclude that only for these orbits 
the use of spectroscopic factors is sensible and unambigu- 
ous. This is in complete accord with the notion that the 
Landau quasiparticle (hole) concept is only valid in the 
immediate vicinity of the Fermi energy [4J-|46| which is 
discussed for nuclei e.g. in Refs. (47|,|4j|. As discussed in 
Ref. [lf|, these quasihole spectroscopic factors are con- 
sistent with the analysis of the (e,e'p) reaction on this 
nucleus for the d 3 / 2 orbit [HI, [5(1]. These references quote 
a result for the lsi/ 2 orbit that is too small and has been 
superseded by more recent experiments [51| . For more 
deeply-bound orbits, including the d 5 / 2l it is much more 
appropriate to consider the complete spectral functions 
as shown in Fig. [2j We note that the very sharp peaks in 
the immediate vicinity of the Fermi energy can also be 
replaced by delta-functions with the strength given by 
the spectroscopic factor. 

By integrating the imaginary part of the propagator 
given in Eq. ([I} for each £j combination up to the Fermi 
energy it is possible to obtain the one-body density ma- 
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FIG. 3. (Color online) Experimental charge density of 
40 Ca Q (solid) compared with the DOM result (dashed). 



trix element 

n tj {r',r) = - [ * dE lmG i:j (r,r';E) 

= (*o|^iOr«l*0>- (20) 

For protons, the point charge distribution is thus ob- 
tained from the diagonal matrix elements of the one-body 
density matrix 



P P ( r ) 



e 

4-7T 



£)(2i + l)n«(r,r). 



(21) 



For a comparison with the experimental charge density 
of 40 Ca it is necessary to fold this distribution with the 
proton charge density. We used the procedure outlined in 
Ref. [52l | which employs 3 gaussians for the proton. The 
mean square radius of the resulting charge distribution 
is obtained from 



<r 2 > 



1 

~Ze 



dr r 2 p ch (r) 



(22) 



and has been employed to constrain the nonlocal HF po- 
tential to generate good agreement with the experimen- 
tal mean square radius of the Ca charge distribution. 
The parameters in Table Q] generate a value of 3.45 fm 
compared to the experimental result of 3.45 fm taken 
from the Fourier-Bessel analysis given in Ref. (42[. We 
compare the calculated charge density with the exper- 
imental one in Fig. [3] It is obvious that there is still 
a significant discrepancy with the experimental charge 
density near the origin which requires further analysis. 
Before addressing this issue in more detail we remark on 



TABLE IV. Radii for proton orbits in 40 Ca for the local and 
nonlocal DOM implementation. 



orbit 



local [fn 



nonlocal [fm] 



Osi/2 
OP3/2 
OP1/2 
0d 5 / 2 

lSl/2 



3/2 



2.34 
2.99 
2.98 
3.54 
3.87 
3.71 



2.36 
2.92 
2.90 
3.36 
3.60 
3.52 



the individual radii of the quasihole orbits and compare 
these with results obtained from the (e,e'p) reaction. By 
constraining the nonlocal HF potential to reproduce the 
mean square charge radius, we obtain radii of individual 
orbits that are somewhat smaller than for the local DOM 
fit. We list the radii for the solutions of Eq. ([6]) (without 
the imaginary part) in Table |IV] and compare the DOM 
results with each other. We first note that the radii listed 
are for point particles and become larger when folded 
with the proton charge density to generate the result for 
the nonlocal case that agrees with the experimental mean 
square radius of the charge distribution. This constraint 
makes the nonlocal radii smaller than the ones for the 
local DOM. The radii for the local DOM are in good 
agreement with the (e,e'p) analysis of Refs. [SiJ [5fl| for 
the valence hole states. A direct comparison of the point 
proton radii is appropriate, since the analysis of the data 
employs bound-state wave functions for point nucleons 
but includes the coupling to an extended charge distri- 
bution of the proton by employing the off-shell (e, p) cross 
section. We do note that the (e,e'p) analysis generates a 
sp wave function with respect to the A — 1 system. This 
is not the case for the calculation of the DOM propa- 
gator with the nonlocal HF potential. Since the DOM 
calculation generates proton scattering wave functions in 
a wide energy domain and also calculates the proton over- 
lap function for removal to valence hole states, it provides 
all the ingredients that are employed in the analysis of 
the (e,e'p) cross sections. Such cross sections, in turn, 
can therefore be used in future DOM calculations to con- 
strain the nonlocal potential to check whether the radii 
obtained for the nonlocal DOM describe the data. In- 
deed, we note that elastic-nucleon-scattering data only 
determine the phase shifts associated with the asymp- 
totic scattering wave functions and do not provide strong 
constraints on the interior scattering waves. 

The possibility of generating the one-body density ma- 
trix from a nonlocal DOM calculation provides access to 
the natural orbits of the system. We discuss these orbits 
before we address the discrepancy of the DOM charge 
density with the experimental one, since the shape of the 
natural orbits provides some initial clues as to what is 
missing in the DOM potentials as implemented so far. 
By diagonalizing the one-body density matrix given in 
Eq. (|20|) one obtains the natural orbits for each £j combi- 
nation together with the corresponding occupation num- 
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r[fm] 

FIG. 4. (Color online) Comparison of the wave functions for 
the d 3 /2 natural orbit (dashed) and the corresponding quasi- 
hole result (solid). The quasihole wave function has a slightly 
larger radius. 

bers. It is therefore possible to write 

ntj(r, r') = MO'), (23) 

n 

with n"Jy , ipntj ( r ) the corresponding occupation numbers 
and wave functions for the natural orbits. We note that 
these wave functions are normalized to unity. While there 
is a correspondence for n-values that are nominally oc- 
cupied between the natural orbits with large occupation 
numbers and overlap functions that correspond to mostly 
occupied states, this is lost for the natural orbits with 
small occupation numbers [Hj]. Indeed, natural orbits 
with small occupation numbers are more confined than 
those with large occupation and have no relation with 
mostly empty orbits associated with, for example, DOM 
potentials. 

We compare in Fig. 0] the proton d 3 / 2 natural orbit 
(solid) with the corresponding quasihole wave function 
(dashed) obtained by diagonalizing Eq. (J6j) at the cor- 
rect energy. Both wave functions are normalized to 1 
and are basically indistinguishable. This feature was also 
observed in the microscopic calculation of Ref. [13 for 
16 0. For the s\/2 orbit two natural orbits are generated 
with large occupation numbers. In this case the com- 
parison in Figs. [5] and [5] shows quite similar behavior. 
The DOM quasihole wave functions exhibit the expected 
shell-model wave function characteristics since they have 
been generated by solving Eq. ^ and occur at quite dif- 
ferent energies. The natural orbits are generated quite 
differently as they require an integration of the spectral 
density over all energies up to the Fermi energy with a 
subsequent diagonalization of the the one-body density 
matrix without any direct reference to a Schrodinger-like 
equation. Nevertheless, the wave functions of natural 



orbits with large occupation numbers appear almost in- 
distinguishable from their quasiparticle counterparts. 

Occupation numbers for natural orbits are collected in 
Table [V] The number of large eigenvalues (comparable 
to unity) corresponds exactly to the expected number of 
the simple shell model. All other eigenvalues are small 
and are associated with wave functions with increasing 
number of nodes. We also include the sum of the oc- 
cupation numbers for each Ij combination in Table [V] 
Somewhat surprising is that the largest deviation occurs 
for the S1/2 orbit (in that case from 2, since there are 
nominally two levels occupied). As observed in Ref. [53j j . 
the largest eigenvalues for the nuclear natural orbits are 
substantially larger than the corresponding one for drops 
of a finite number of 3 He atoms [54j ■ Since short-range 
correlations associated with the underlying bare inter- 
action are included in the work of Refs. [53[ and [54j . 
this difference is mostly related to the much stronger re- 
pulsion between 3 He atoms which e.g. in the liquid at 
saturation leads to a depletion of the Fermi sea of more 
than 50% j55|. Nucleon-nucleon interactions typically 
generate 10-15% depletion due to short-range and tensor 
correlations [56| . The inclusion of short-range correla- 
tions is partly accomplished by the assumed energy de- 
pendence of the volume term of the imaginary part of the 
DOM potential which is based on nuclear matter calcula- 
tions 57]. Such an imaginary term is responsible for the 
global depletion of orbits including those that are deeply 
bound. We note that the complementary admixture of 
high-momentum components is not yet incorporated by 
the current DOM implementation, as will become clear in 
the following discussion. We also observe that the large 
occupation numbers of the natural orbits calculated for 
16 in Ref. H are 5-10% larger than the DOM num- 
bers quoted in Table [V] It appears reasonable to inter- 
pret this difference to be due to the proper inclusion of 
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r[fm] 



FIG. 5. (Color online) As in Fig. [4] but for proton S1/2 orbit 
with one node. 



10 



TABLE V. Occupation numbers of natural orbits. 



n 


Sl/2 


P3/2 


Pl/2 


^5/2 


d.3/2 


h/2 


fs/2 


1 


0.926 


0.921 


0.905 


0.899 


0.858 


0.109 


0.064 


2 


0.881 


0.072 


0.062 


0.037 


0.032 


0.024 


0.020 


3 


0.032 


0.021 


0.020 


0.015 


0.014 


0.010 


0.010 


4 


0.015 


0.010 


0.009 


0.007 


0.007 


0.006 


0.005 


5 


0.007 


0.005 


0.005 


0.004 


0.004 


0.003 


0.003 


y 


1.86 


1.03 


1.00 


0.96 


0.92 


0.15 


0.10 



t • 1 • r 




r [fm] 



FIG. 6. (Color online) As in Fig. |4]but the proton S1/2 orbit 
with no node. 




k [fm 1 ] 

FIG. 7. (Color online) Comparison of the total momentum 
distribution calculated according to Eq. (|24[) (dashed) with 
the one obtained from the quasihole contributions (solid). 



low-energy, and therefore long-range correlations that are 
represented by the surface components of the DOM po- 
tentials. The imaginary part of the DOM potential up 
to about 50 MeV is dominated by the surface potential. 
The associated absorption is well-constrained by the dif- 
ferential cross sections for elastic nucleon scattering in 
this energy domain. It is then also unavoidable that sp 
strength from below the Fermi energy is removed to this 
energy domain when such potentials are employed in the 
solution of the Dyson equation. 

DOM implementations usually generate occupation 
numbers for quasihole and other bound orbits Q based 
on approximate expressions. It is therefore useful to com- 
pare these with the more accurate results obtained by 
integrating the corresponding strength up to the Fermi 
energy. This can also be accomplished by using the quasi- 
hole wave functions and performing an integration over 
r and r' involving the one-body density matrix. For the 
proton <i 3 /2 quasihole orbit we find an occupation number 
of 0.86 which is identical to the one for the natural orbit. 
This is hardly surprising since these wave functions are 
almost identical as shown in Fig. @] The DOM result of 
Ref. [18| yields 0.82 reasonably close to the Green's func- 
tion result of the nonlocal DOM. For the lsi/2 quasihole 



state calculated at -9.8 MeV we obtain an occupation 
number of 0.88 a little larger than the result obtained 
with the local DOM which generates 0.85. For the 0si/2 
state the corresponding numbers are 0.93 (nonlocal) and 
0.93 (local) in complete agreement. This suggests that 
the corresponding approximate expression in the local 
DOM may be appropriate for the occupation numbers 
deep in the Fermi sea. It is however well known that 
for such orbits unreasonable spectroscopic factors can be 
generated. For the Osi/2 orbit the local DOM yields for 
example an unphysical spectroscopic factor of 1.11 (see 
Table IIII[) . The corresponding calculation for the non- 
local DOM yields 0.98, a not unphysical result but also 
not useful since it is larger than the properly calculated 
occupation number of 0.93. This is hardly surprising 
since the calculation employing Eq. © neglects the role 
of the imaginary part of the self-energy completely. As 
discussed earlier, it is therefore only useful to consider 
spectroscopic factors near the Fermi energy where the 
imaginary part of the self-energy is insignificant or zero. 
At other energies, it is more appropriate to consider the 
spectral function as illustrated in Fig. [2j 

As the discussion of natural orbits has shown, the effect 
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of short-range and tensor correlations are only included 
in so far as orbits below the Fermi energy are depleted 
by the presence of a substantial imaginary part at large 
positive energy associated with the volume contribution 
to the DOM self-energy. The associated presence of high- 
momentum components is not explicitly guaranteed and 
we will analyze this in the following. The total proton 
momentum distribution for protons (normalized by Z) is 
obtained from 



(24) 



We obtain the partial momentum distributions nj>j(k) by 
first generating the momentum-space spectral function 
by performing a double Fourier-Bessel transform of the 
spectral density in coordinate space 



oo />oo 



JO 



S ej (k;E) = / drr 2 l dr 1 V 2 jt(kr)lmG ej (r,r' '■ E)j t (kr'). 

The momentum distribution for a given £j is then ob- 
tained from 



n tJ (k)= / dES tj (k;E). 



(26) 



In Fig. [7] we display the total proton momentum distribu- 
tion by the dashed line. For comparison we also show the 
momentum distribution from the quasihole wave func- 
tions (normalized to one) by the solid line. As dis- 
cussed in previous work (see e.g. Refs. [H], |3f| l53|). 
these quasihole contributions are mostly associated with 
wave functions near the Fermi energy and hardly contain 
any high-momentum components. The presence of high- 
momentum components is demonstrated by the dashed 
line in Fig. [7J We emphasize their contribution by show- 
ing in Fig. [8] the momentum distribution weighted by 
k 2 . Interestingly we find about 10% of the protons ac- 
tually have momenta beyond 1.4 fm~ . This number is 
in reasonable agreement with the 10 % g enerated for 16 
in the calculations of Refs. [13, HIJ l53| These calcula- 
tions generate high-momentum components that are in 
quite good agreement in the aggregate with the results of 



Ref. 33] . It is therefore clear that the present version of 



the DOM includes sufficient flexibility to represent these 
experimentally well established ingredients at least in ag- 
gregate. 

Looking in more detail we note that the expected be- 
havior of high- momentum components, i.e. increasing 
importance with increasing separation energy, is not con- 
tained in the DOM spectral functions. We illustrate this 
observation in Fig. |9]by plotting the d^/2 spectral func- 
tion in momentum space at different energies; starting at 
-25 MeV with steps of 25 MeV all the way to -150 MeV. 
The figure illustrates that the shape of the momentum 
content of the spectral function hardly changes as a func- 
tion of energy, especially when momenta above 1.4 fm 
are considered. This latter feature is completely oppo- 
site to the effect expected of short-range correlations. As 
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FIG. 8. (Color online) Momentum distribution for protons as 
in Fig. [7] but weighted by k 2 . 



discussed at length in Refs. [3J, |35|, (53|, the presence of 
high-momentum components becomes more pronounced 
with decreasing energy (away from the Fermi energy) un- 
like the results shown in Fig. [9l The former result can 
be easily understood on the basis of simple considerations 
involving momentum conservation and the location of the 
relevant 2hlp states that are required for the admixture 
of high- momenta [3(| ■ We note that this effect has been 
experimentally confirmed by the absence of appreciable 
high-momentum components in the valence hole states in 
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k [fm ] 



FIG. 9. (Color online) Momentum-space spectral function 
for d 3 / 2 quantum numbers at different energies. The highest 
curve is obtained at -25 MeV and each successive lower curve 
(at small momenta) represents a 25 MeV step lower in energy 
with the last curve representing the spectral function at -150 
MeV. 
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FIG. 10. (Color online) Momentum-space spectral function 
for S]/ 2 quantum numbers at different energies as explained 
in the text. 



208 Pb [58j]. The results of Ref. j33| further illustrate that 
high-momentum components emerge with decreasing en- 
ergy and dominate at energies substantially below the 
bottom of the traditional potential well used to describe 
mean-field nucleons. In Fig. [10] we display the momen- 
tum content for the S\/2 orbit at the same energy val- 
ues. In this case the short-dash-dot curve at -25 MeV is 
reminiscent of the lsi/2 quasihole wave function and the 
long-dash-dot curve at -50 MeV is close to the quasihole 
peak of the Osi/2 orbit. At lower energy this shape per- 
sists but the high-momentum content (apart from slowly 
decreasing) exhibits no essential change in energy as for 
the d 3 /2 channel. 

In order to describe the correct behavior of the high- 
momentum components in the DOM it will be necessary 
in the future to make the geometry of the potential de- 
pendent on energy. Indeed, by reducing the radius of 
the confining nuclear potential with decreasing energy, 
one may expect to raise the high-momentum content and 
generate the behavior predicted in Refs. [13, and ex- 
perimentally confirmed in (33j . Since the geometry of 
the DOM potential has been assumed independent of en- 
ergy in the current implementations, this will increase 
the computational effort substantially since the applica- 
tion of the subtracted dispersion relation will have to be 
performed also as a function of the coordinates for which 
the real par t of the dispersive part is required. The work 
of Refs. |34l . Hfl was performed in momentum space and 
it may be necessary to consider DOM implementations 
which rely on momentum-space formulations, at least as 
far as short-range correlations are concerned. 

Having established some missing ingredients in the de- 
scription of high-momentum components, we now argue 
that this has consequences for the description of the nu- 
clear charge density. As discussed in Ref. |59j, the role 



of short-range correlations is to remove some nuclear 
charge, present in the mean-field description in terms of 
the occupied s^n states, from the origin to larger radii 
but not to the surface, which is dominated by quasihole 
contributions. While some of this charge returns to the 
origin as partially occupied higher S\/2 states, most of 
this strength is associated with higher ^-values, similar 
to the results obtained in Refs. [35j, |53|. It is therefore 
reasonable to expect that a proper treatment of short- 
range correlations with the attendant presence of high- 
momentum (higher £) components (constrained by the 
experimental data 33]) will make it possible to obtain 
an accurate fit to the nuclear charge density in a DOM 
framework. 

It is well known that the sp propagator allows for the 
calculation of the energy per particle from the contri- 
bution of the underlying two-body interaction. For the 
present case, it is useful to employ this result in momen- 
tum space. The energy per proton of the ground state 
can, for example, be obtained by calculating [27| 

£( 40 Ca) 1 . f°° lll2 k 2 „. 

i roo pep 

+^B 2 J + 1 ) / dkk2 / dE ES ej (k;E). (27) 



For the present DOM potential we obtain only -2.91 MeV 
per proton which includes the effect of the Coulomb in- 
teraction. A similar calculation for the neutrons yields 
-6.51 MeV per neutron for a total of -4.71 MeV per par- 
ticle. This result represents about 60% of the experimen- 
tal result. This is a remarkable result since the spectral 
information and the location of the bound levels in com- 
bination with a considerable wealth of elastic scattering 
data is described by the DOM self-energy. However, also 
in this case we can point to the lack of the correct descrip- 
tion of high-momentum components that can resolve this 
issue. In Ref. 35] it was shown that the quasihole con- 
tribution to the energy per particle is about 35% in 16 
whereas 65% is generated by the continuum contribution 
at large negative energies where high-momenta dominate. 
This result is noteworthy also since only 10% of the nu- 
cleons are considered to have high momenta as confirmed 
by experiment. A similar situation appears to apply in 
the case of the DOM analysis of 40 Ca. Since the total 
number of high-momentum components appears reason- 
able, it appears that their appearance at more negative 
energy will be able to resolve part of the discrepancy for 
the total energy of the ground state. It must also be 
noted that an important contribution from three-body 
forces may have to be considered. It appears therefore 
reasonable to expect that all data that are not yet well 
reproduced at present, can be better described in a future 
DOM implementation which incorporates the contribu- 
tion of about 10% of high-momentum nucleons with the 
correct energy dependence. 
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IV. SUMMARY AND CONCLUSIONS 

The present work aims at extending the DOM ap- 
proach, which so far has been mostly applied to describe 
elastic nucleon scattering, into the domain below the 
Fermi energy by employing additional experimental data 
to constrain the potentials. By introducing an explicit 
nonlocal HF-like potential, it is possible to reinterpret 
the DOM potential as a nucleon self-energy when the 
corresponding nonlocality correction is implemented to 
generate the intended normalization of the DOM poten- 
tial containing the local but energy-dependent HF con- 
tribution. This procedure has been adopted for DOM 
potentials that were previously obtained for 40 Ca [l8j]. 
The possibility to interpret the DOM potential below 
the Fermi energy as the nucleon self-energy broadens the 
links with experimental data substantially. The solution 
of the Dyson equation below the Fermi energy with this 
self-energy then leads to the nucleon sp propagator and 
the corresponding one-body density matrix. A Perey- 
Buck type nonlocal potential was chosen to represent the 
HF potential. Its parameters were chosen to describe the 
energies of the valence hole states and the mean square 
radius of the charge-density distribution. The latter fea- 
ture illustrates the new possibility to constrain the DOM 
potential by data that pertain to information associated 
with one-body properties of the nuclear ground state, 
since the sp propagator provides this access for any one- 
body operator. Various quantities that are obtained from 
approximate expressions with the usual local form of the 
DOM potential are compared with the nonlocal solution. 
Spectroscopic factors near the Fermi energy appear to 
be stable quantities, but are no longer useful for deeply 
bound states. Instead, we advocate the construction of 
the complete spectral function in particular when com- 
parison with nucleon knockout experiments are consid- 
ered. The nonlocal HF potential also limits the binding of 
the lowest Sx/a orbital in agreement with corresponding 
experimental information. This orbit tends to become 
too deeply bound with the local version of the DOM po- 
tential. 

The diagonalization of the one-body density matrix al- 
lows for the study of natural orbits and associated occu- 
pation numbers (eigenvalues). Results are qualitatively 
similar to microscopic calculations performed earlier for 
16 and drops of a finite number of 3 He atoms. For 
each orbit that is filled in a simple mean-field picture, 
there is a corresponding large eigenvalue, while all other 



eigenvalues for this £j combination are at least an or- 
der of magnitude smaller. Independent of whether one 
or two such orbits are occupied below the Fermi energy, 
like for the and Sx/2, respectively, there is little dif- 
ference between the quasihole (overlap) and natural orbit 
wave functions. Although the mean square radius of the 
charge distribution agrees with the experimental value, a 
comparison with the complete density distribution shows 
that too much charge is calculated near the origin. By 
studying the momentum content of the spectral func- 
tion and associated momentum distribution, we observe 
that current DOM potentials generate about 10% high- 
momentum com pon ents in agreement with experimen- 
tal observations [33[ for light nuclei. Microscopic calcu- 
lations of short-range correlations however demonstrate 
that high-momentum components become increasingly 
important with decreasing energy away from the Fermi 
energy, which is also confirmed by experiment. This fea- 
ture must be represented by a change in geometry of the 
potential with decreasing energy, which is currently not 
present in DOM potentials but can be implemented in fu- 
ture applications. Since short-range correlations remove 
charge from the origin (sx/2 orbits) and place it at larger 
radii although not at the surface [59| , it is expected that 
the inclusion of high- momentum components in the DOM 
potentials will allow a better description of the complete 
charge distribution. Similar considerations suggest that 
these high-momentum components will also generate a 
substantially improved energy per particle. 

Since the DOM generates both scattering wave func- 
tions at positive energy as well as quasihole overlap func- 
tion, it will in the future also be possible to describe 
(e,e'p) cross sections directly and use these data to con- 
strain the DOM potential further. At present, the indi- 
rect comparison with spectroscopic factors derived from 
these data suggests no inconsistency. Nevertheless, it is 
important to explore the interior scattering wave func- 
tion of the outgoing proton in this process, since it may 
be sensitive to the nonlocal features of the potential. 
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